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Abstract 

The behaviour of baby universes has been an important ingredient in understand- 
ing and quantifying non-critical string theory or, equivalently, models of two- 
dimensional Euclidean quantum gravity coupled to matter. Within a regularized 
description based on dynamical triangulations, we amend an earlier conjecture 
by Jain and Mathur on the scaling behaviour of genus-^f surfaces containing par- 
ticular baby universe 'necks', and perform a nontrivial numerical check on our 
improved conjecture. 
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1 Introduction 



In physics, the intriguing concept of baby universes as the offspring of a parent 
universe - acquiring a hfe of their own while remaining connected to the latter 
- has been most concretely realized in 2d Euclidean quantum gravity models. 
In these models the partition function for a compact genus-gi surface coupled to 
conformal matter with central charge c is given by 

=^T(^)-V^, (1) 

a formula which is exact at the conformal point (see e.g. [1]). In (1), A denotes 
the area of the surface, /i a bare (cutoff-dependent) cosmological constant, and 
7(51) the so-called susceptibility exponent or entropy exponent, given by 

7(^) = 7o + ^(2-7o), 7o-7(0)- '"'"^^'^~'^^'"''^ - (2) 

For c = 0, where no matter is present, the partition function (1) has an entropic 
interpretation, namely, it "counts" the number of surfaces of genus g with area 
A. After introducing a regularization of the continuum theory, for example, by 
using dynamical triangulations (DT), this counting becomes literal. The con- 
tinuum area is then A = \/3a^N/A, where is the number of triangles, and a 
the edge length of the equilateral triangles used in DT, simultaneously acting 
as a UV-cutoff of the theory. The number of distinct piecewise linear surfaces 
(triangulations) of genus g one can construct from N such triangles is 

Z^a\N) = k{g) Ar7(5)-3e/^oiv ^ 0{1/N)) , (3) 

where /Iq is independent of the genus g and k{g) is a constant. 

In an influential paper, Jain and Mathur (JM) related the entropy exponent 

7o to the creation of baby universes [2]. They showed that the average number of 
so-called "minbus" (minimal neck baby universes of disc topology) of regularized 
area B is given by 

nA{B) = k{g) A{1- B/Ay^'^'^ B^°-\ (4) 

where A is the regularized area of the total surface, and we assume B < A/2. 
Cutting the surface open along the minimal neck will produce two disconnected 
surfaces, a genus-si surface of area A — B with one (S^-boundary (the minimal 
neck), and a baby universe of area 5, whose topology is a sphere with a single S*^- 
boundary (the same minimal neck), that is, a disc. In the following we will always 
work in a suitable DT-ensemble of triangulated surfaces and use the notation A 
and B for the number of triangles, thus ignoring the factor a^-\/3/4 which relates 
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the number of triangles to the area of the surface. Likewise in discrete units, the 
length of a minimal neck can be 1, 2 or 3 links (= triangle edges), depending on 
the choice of DT-ensemble, i.e. the regularity conditions imposed on the triangle 
gluings. 

The leading power-law behaviour of B in (4) gives us a convenient way to 
determine 70 in numerical simulations and has been used extensively in two di- 
mensions^, beginning with [3], as well as in higher-dimensional gravity, starting 
with [4]. 

While the derivation of (4) is robust and correct, JM put forward generalized 
relations, based on a certain conjectured formula involving non-minimal necks 
(eq. (9) in [2]), which also divide a surface into two pieces A — B and but have 
an arbitrary length £ < \/~A. The main purpose of this letter is to demonstrate 
that the formula conjectured in [2] is incorrect. We will discuss how to modify it 
appropriately, while retaining many of the results derived there. Our key obser- 
vation will be that the necks of baby universes are by construction rather special 
curves on the two-surface, whose scaling behaviour (as function of the surface 
area A) is different from that of generic curves. We will also perform a nontrivial 
numerical check of our improved conjecture and suggest further applications of 
our new result. 

2 The conjecture and how to modify it 

In this section we will limit ourselves to the case c = of 2d gravity without 
matter, leaving the discussion of general central charge c < 1 to Sec. 4. 

Consider a (triangulated) surface of area A and spherical topology, except for 
n boundary loops of length £j, i = 1, . . . , n, counting the number of edges in the 
boundaries. Assume also that none of the boundary loops can be deformed into a 
loop of shorter length in the same homotopy class, unless the deformation sweeps 
an area which is a sizeable fraction of A. JM conjectured that the number of such 
surfaces behaves like 

n 

ZMA) ~ ^(^=°)(A) A" J] £7(1+^"). (5) 

j=i 

The factor A" Z^^^'^^i^A) is uncontroversial, with A"' counting the number of ways 
n boundary loops can be located on a surface of area A. In addition, JM made the 
general ansatz where i is less than some fraction of \fA. To determine 

^ "Integrating out" baby universes in two dimensions [5] leads to an alternative theory of 
2d quantum gravity known as "causal dynamical triangulations" (CDT) [6]. By introducing a 
coupling constant for baby universe creation one can relate the two models [7]. 
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a, they calculated the genus-one partition function from Z2{A,i,i) by gluing 
together the two boundaries to form a torus and integrating over £ up to \/A, 




dii Z2{A,i,i), 



(6) 



where the factor £ in the integrand takes into account the number of different 
gluings of the two boundary loops. Comparing to (1) and (2) one obtains a = 70. 

A number of interesting results were derived in [2] using the ansatz (5). For 
example, it is straightforward to obtain the relation '~f{g) = 70 + (?(2 — 70) by 
employing (5) for a surface A with g pairs of boundaries of equal loop length ii, 
i — 1, . . . ,g, and thus n — 2g. General genus-g' surfaces are created by gluing 
together the loops of each pair and integrating over the corresponding length 
variable ii, generalizing (6). It is easy to see that each such operation will create 
a factor A'^^'^o. 

As pointed out in [2], the number of spherical surfaces of area A with n 
boundaries which cannot be deformed without increasing their length is strictly 
smaller than the number of surfaces when the n boundaries are arbitrary, for 
which standard formulas exist^. This happens because boundaries which serve 
as baby universe necks in the sense of JM are special curves. However, the con- 
sequences of this are even more drastic than envisaged in [2], and are directly 
related to why the ansatz (5) is not entirely correct. The boundary loop in (5) is 
by definition a geodesic curve, and it is well known that geodesies scale anoma- 
lously in the DT ensemble of surfaces [9, 10, 11, 12, 13]. In fact, the dimension 
of geodesic curves is volume^^^ and not volumc^/^, as one might have expected 
naively, implying that the HausdorfF dimension dh of surfaces in the DT ensemble 
is 4. This is reflected in the scahng behaviour of the expectation values (w.r.t. 
the ensemble average) 



{R)a--A^/'^\ (A(r))A~r'^N r</l^/'^^ where 4 = 4 (c = 0), (7) 



where R is the linear extension of a surface of area A, and A[r) the area contained 
within a geodesic distance r from a given point. 

The necks of length £ along which JM cut surfaces into disconnected pieces are 
geodesic curves, which means that their ensemble average is much shorter than 

^What could have alerted the authors of [2] to the fact that they were not quite on the right 

track is that their formula for n > 3 and c = is identical to the standard expression (without 
any restrictions on the loops) derived in [8], namely, 




This formula is valid also for ii > y/A and shows that the cutoff in £ for ordinary boundaries is 
of order y/A, as one would expect from the dimensionalities of £ and A. 
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the generic ^J~A used in [2] as upper limit in integrations like (5). (Note that the 
main contribution to this integral comes precisely from the upper limit.) Instead, 
according to (7), A^l'^ should be used as upper limit, and, more generally, A}!'^^ if 
the Hausdorff dimension is dh- An alternative ansatz, which will reproduce most 
of the results of [2] , is to replace (5) by 

n 
i=l 

This formula is supposed to be valid for < A^^'^f' , with the understanding that 
any integrations over £j arc to be performed from the minimal neck length up to 
^i/dh por £i > A^^'^^ the function Z'^{A,£i) is assumed to vanish fast. 



3 Testing the new ansatz 

As an application of the new ansatz (8), let us reconsider the torus case, where 
relation (6) is now replaced by 

We have performed a nontrivial quantitative test of the ansatz (8) by measuring 
the probability distribution Pa{^) of the length t of the shortest non-contractible 
loop in the DT ensemble of triangulations of the torus with A triangles. According 
to (8), the distribution should be 



tZ'^{A,£,tj I 



where F(0) = 1 and F{x) goes rapidly to zero when x> \. 

In practice, to find shortest non-contractible loops we need a method to de- 
termine whether a given loop is contractible or not. By a different method, we 
first constructed two loops ai, 0:2 generating the fundamental group of the torus. 
Whether an arbitrary loop is contractible can then be established by measuring 
its intersection numbers with these two generators. It is contractible if and only 
if both intersection numbers vanish. The shortest non-contractible loop based at 
any given vertex v can be found by performing a so-called breadth-first search 
along the edges starting at v. Encountering a vertex that has been visited before 
means that one has found a loop in the triangulation. The first such loop one 
meets which is non-contractible will automatically have minimal length. Since 
any non-contractible loop must intersect at least one of the ctj, in order to find 
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Figure 1: Plot of a rescaled version A'^/^PAii)/{e/A'^/^) of the probability distri- 
bution of the length of the shortest non-contractible loop on the DT ensemble of 
triangulated tori, as function of the rescaled length x = i/A^^^, with £ G [2,21] 
and areas A in the range between 125 and 64000. 



the overall shortest non-contractible loop it suffices to repeat the above procedure 
for all vertices contained in either of the generating loops aj.^ 

Fig. 1 shows the superimposed measurements of the rescaled probability dis- 
tribution y4^/^P^(^)/(£/A^/'^) for various values of i and A. The computer simu- 
lations were performed in the DT ensemble consisting of triangulations dual to 
connected 0^-graphs with A vertices and toroidal topology, including self-energy 
and tadpole graphs. (To be precise, the simplest self-energy subgraphs consisting 
of two lines and tadpole subgraphs consisting of a single line were left out to 
ensure that all vertices in the triangulation have an order larger than or equal to 
three. This is not expected to affect the final result in any way.) 

If the ansatz (8) is correct we should observe a universal function of a; = i/A^^'^ 
which goes to a constant for a; — )■ 0. Inspection of Fig. 1 reveals that the expected 
finite-size scaling is indeed satisfied with good accuracy, and we obtain a universal 
function F{x). By contrast, the computer simulations are in clear disagreement 

^To optimize performance one should choose generating loops which are relatively short. In 
line with the result presented in this paper, this means that their length is of the order V^^'^'^ , 
implying an expected run-time of the full algorithm of the order — y^-^^. 
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with the ansatz of [2] which predicts a (^-independent) probabihty distribution 

py^)(£)^_^ for l<^^A, (11) 

as noted in [14]. Our improved ansatz also resolves a discrepancy found in [14], 
namely, that (11) leads to an average loop length {t) ~ A^l'^ ^ while one would ex- 
pect (£) ~ A^l'^ because of the anomalous scaling dimension of geodesic distance. 
It is exactly the latter relation which follows from our suggested probability dis- 
tribution (10), based on the ansatz (8). 



4 2d gravity coupled to conformal matter 

The arguments presented above for c = need to be refined for the general case 
c < 1, when matter is coupled to 2d gravity. The idea put forward in [2] of 
how to tackle this situation is to integrate out all matter degrees of freedom 
and work with the resulting partition function (which looks purely geometric, 
like the one for c = 0), and then simply repeat the cutting-and-gluing arguments 
for baby universes used in the matter-free case. Although this is a well-defined 
procedure for a single closed surface of area and will lead to a partition function 
Z(f)(A), this is no longer true in the presence of boundaries where, to start with, 
boundary conditions for the matter must be specified. Examples are so-called 
"free boundary conditions" (where the matter configurations on the boundary are 
integrated over), or other specific prescriptions like Dirichlet boundary conditions 
for the case of a scalar field, say. No matter what choice one makes, additional 
weights associated with the matter on the boundary will arise, which are simply 
not taken into account in the treatment of [2], which crucially features "baby 
universe surgery" along such boundaries. 

Let us illustrate this point by the case of Ising spins coupled to DT, with spins 
located at the centres of triangles. The Ising model has a critical point when 
defined on the DT ensemble, describing a c = 1/2 conformal field theory coupled 
to 2d gravity [15]. If we consider a surface containing a neck and an associated 
baby universe, there is clearly an energy associated with the spin configurations 
of neighbouring triangles across the loop forming the neck, which will appear as 
a weight in the statistical sum over spin configurations. There is no obvious way 
to recover this energy contribution exactly from the separate partition functions 
of the two disjoint surfaces after cutting open the initial surface along the neck. 

A perhaps more pertinent question is when and whether such unaccounted 
boundary contributions to the energy will matter. In the case of minbus the 
boundaries in question will have a minimal number of edges (one, two, or three, 
depending on the choice of DT ensemble) , and any energy contribution associated 
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with them will become negligible in the continuum limit. This explains why the 
relations derived by JM are robust for the case of baby universes of minimal neck 
size, even when c 7^ 0. As already mentioned in passing, their relation (4) for the 
average minbu number has been used to measure 70 also in the matter-coupled 
case, and agreement with the theoretical value for 70 was found. 

The derivations in [2] building on the ansatz (5) for general, non-minimal 
necks are much harder to justify. When the boundary lengths I become of the 
order of the linear system size (it is precisely those large ^-values which give 
important contributions to integrals like (6), as noted earlier), the contributions 
from boundary energies become potentially significant, but are very difficult to 
control. Nevertheless, let us for the sake of simplicity assume that we can ignore 
such boundary energies when decomposing surfaces and that a suitable general- 
ization of the ansatz (8) is valid, and see where this leads us. 

For definiteness, consider the Ising model on the DT ensemble as described 
above. As long as the coupling constant of the Ising model stays away from its 
critical value, the geometric behaviour of the model will be identical to that of 
the pure-gravity case with c = and = 4, and captured by our earlier ansatz 
(8), together with the new insight that the linear size of the surface is not given 
by the 'naive' A^^'^, but by A^^^. By contrast, if we fine-tune the Ising couphng 
to its critical value, the interaction between geometry and spins will be such 
that the fractal structure of spacctime is altered, manifesting itself in a larger 
HausdorfF dimension dh{Ising) >4. This also implies that the effective linear size 
of the system is decreased to A^/'^h{ismg) ^ makes it tempting to conjecture that 
a straightforward generalization of (8), namely, 

n 

i=l 

is a suitable ansatz in the matter-coupled case too. It accommodates the change 
in Hausdorff dimension, as well as reproducing the correct partition function 
2'(f=^)(/l) for the Ising model on a torus. The general formula for as a function 
of the central charge c of the conformal field theory to be used in (12) is [16] 



4(c) = 2^^t^, (13) 

together with the expression given in (2) for 70(c). The relation which generalizes 
(10) is 

Ml i ^ \\lo{.c)\dh{c)-l ( t \ 



where -F(O) = 1 and Fix) falls off fast for a; > 1 (and recall that 7o(c) < 0). We 
conclude that A^I'^^'^'^^PjiXt) should be a universal function x°^F{x) of the rescaled 
length variable ij A^l'^^^"^ , with a = |7o(c)|dft(c) - 1. 
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One could test the prediction (12) numerically for the Ising model, where a 
should be 0.40. A more clear-cut test may be topological 2d gravity, correspond- 
ing to c=— 2, for which 70 = — 1, (i/i = 3.56 and thus a = 2.56. A special feature of 
this case is that triangulated surfaces of spherical topology can be constructed re- 
cursively with the correct weight, removing the need for Monte Carlo simulations 
[17]. This was employed in [18] to check (13) numerically with great precision. 
If one can generalize the recursive generation of triangulated surfaces to toroidal 
topology it may be possible to obtain a high-precision determination of a too. 

Lastly, there exists a formula for dh{c) different from (13), namely 



dating all the way back to one of the original papers of quantum Liouville theory 
[19]. It agrees with (13) for c = 0, where dh — 4:. However, it gives d/j — >■ 
for c — > — 00 (where one would expect dh — 2), and — > 00 for c — > 1. It has 
the interesting feature that the exponent a introduced above equals 1, not only 
for c = 0, as we have already verified, but for all c < 1. By now it is believed 
that the dh{c) of eq. (13) does not reflect directly any geometric properties of the 
fluctuating geometry, but rather critical aspects of the matter flelds coupled to 2d 
gravity. This was made explicit in [20] , where it was shown that the spin clusters 
of the matter flelds had exactly the dimension (15) for both c=— 2 and c — 1/2, 
while measurements of the geometric d^ clearly differed from (15). Measuring 
the exponent a for c= —2 and c= 1/2 as mentioned above would give us a new 
method, compared to [20], to decide whether (13) or (15) represents the geometric 
Hausdorff dimension of 2d quantum gravity coupled to matter. 
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